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Abstract 

The measurement of the mixing angle #13, sign of Am 2 3 and the CP or T 

X' 

violating phase 5 is fraught with ambiguities in neutrino oscillation. In this 
paper we give an analytic treatment of the paramater degeneracies associated 
with measuring the i/„ — > v e probability and its CP and/or T conjugates. For 
CP violation, we give explicit solutions to allow us to obtain the regions where 
there exist two- fold and four- fold degeneracies. We calculate the fractional 
differences, (AO/ 6), between the allowed solutions which may be used to 
compare with the expected sensitivities of the experiments. For T violation 
we show that there is always a complete degeneracy between solutions with 
positive and negative Am 2 3 which arises due to a symmetry and cannot be 
removed by observing one neutrino oscillation probability and its T conjugate. 
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Thus, there is always a four fold parameter degeneracy apart from exceptional 
points. Explicit solutions are also given and the fractional differences are 
computed. The bi-probability CP/T trajectory diagrams are extensively used 
to illuminate the nature of the degeneracies. 
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I. INTRODUCTION 



The discovery of neutrino oscillation in atmospheric neutrino observation in Super- 
Kamiokande p] and the recent accumulating evidences for solar neutrino oscillations 
naturally suggests neutrino masses and lepton flavor mixing. It is also consistent with the 
result of the first man-made beam long-baseline accelerator experiment K2K [[| . Given the 
new realm of lepton flavor mixing whose door is just opened, it is natural to seek for a pro- 
gram of exploring systematically the whole structure of neutrino masses and lepton flavor 
mixing. 

Most probably, the most difficult task in determining the structure of lepton mixing 
matrix, the Maki-Nakagawa-Sakata (MNS) matrix [|J, is determination of CP violating 
phase 5 and a simultaneous (or preceding) measurement of \U e s\ = sin 6*13. We use in this 
paper the standard notation for the MNS matrix with Am^ = m| — m 2 . The fact that 
the most recent analyses of the solar neutrino data || strongly favor the large mixing angle 
(LMA) Mikheev-Smirnov-Wolfenstein (MSW) solution || is certainly encouraging for any 
attempts to measure leptonic CP violation. 

Since we know that #13 is small, sin 2 2^13 < 0.1, due to the constraint imposed by the 
Chooz reactor experiment |7j and we do not know how small it is, there is two different 
possibilities. Namely, (A) # 13 is determined prior to the experimental search for leptonic 
CP violation, or (B) not. The case (A) is desirable experimentally. To determine unknown 
quantities one by one, if possible, is the most sensible way to proceed with minimal danger of 
picking up fake effects. But since there is no guarantee that the case (A) is the case, we must 
prepare for the case (B). Even in the case (A), experimental determination of 6 I 13 always 
comes with errors, and one must face with the similar problem as in the case B within the 
experimental uncertainties. Moreover, it is known that determination of #13 in low energy 
conventional super-beam type experiments suffers from additional intrinsic uncertainty, the 
one coming from unknown CP violating phase 5. See Ref. || for further explanation and a 
possible way of circumventing the problem. Therefore, the determination of S and #13 are 
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inherently coupled with one another. 

Even more seriously, it was noticed by Burguet-Castell et al. || that there exist two sets of 
degenerate solutions (Si, 9\ 3 ) (i=l, 2) even if oscillation probabilities of P(v^ — > v e ) = P(v) 
and its CP conjugate, CP[P(i/)] = P (v^ — > u e ) is accurately measured. They presented 
an approximate but transparent framework of analyzing the degeneracy problem, which 
we follow in this paper. It was then recognized in Ref. |TU[ that unknown sign of Am^ 
leads to a duplication of the ambiguity, which entails maximal four-fold degeneracy (see 
below). It was noticed by Barger et al. [JO]] that the four- fold degeneracy is further multiplied 
by an ambiguity due to approximate invariance of the oscillation probability under the 
transformation # 2 3 — * f — 6*23- A special feature of the degeneracy problem at the oscillation 



maximum was noted and analyzed to some detail |p|,p] . Recently, the first discussion of the 
problem of parameter degeneracy in T violation measurement is given in Ref. |TI| . 

Meanwhile, there were some technological progresses in analyzing the interplay between 
the genuine CP phase and the matter effects in measuring leptonic CP or T violation in neu- 
trino oscillation, the issue much-discussed but still unsettled |IB,|T1|. The authors of Refs. [10 



and |12| introduced, respectively, the "CP and T trajectory diagrams in bi-probability space" 
for pictorial representation of CP-violating and CP-conserving phase effects as well as the 
matter effect in neutrino oscillations. They showed that when these two types of trajectory 
diagrams are combined it gives a unified graphical representation of neutrino oscillations 
in matter [ 12| . We demonstrate in this paper that they provide a powerful tool for un- 
derstanding and analyzing the problem of parameter degeneracy, as partly exhibited in 



Refs. 10.15.11 



It is the purpose of this paper to give a completely general treatment of the problem of 
parameter degeneracy in neutrino oscillations associated with CP and T violation measure- 
ments. We elucidate the nature of the degeneracy, and determine the region where it occurs, 
namely, the regions in the P-CP[P] (and P-T[P]) bi-probability space in which the same- 
Amf 3 -sign and/or the mixed- Am^-sign degeneracies take place. While partial treatment 
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of the parameter degeneracy has been attempted for CP measurement before ||8|-|TTf such 
general treatment is still lacking. We believe that it is worthwhile to have such an overview 



of the parameter degeneracy issue to uncover ways of resolving this problem. See []16-19 
for recent discussions. 

We present the first systematic discussion of parameter degeneracy in T measurement 
following our previous paper in which we set up the problem [12]. We uncover a new 
feature of the degeneracy in T measurement. Namely, we show that for a given T trajectory 
diagram there always exists an another T diagram which completely degenerates with the 
original one and has opposite sign of Amf 3 . It means that for any given values of P(u) and 
T[P(i/)] there is two degenerate solutions of (5, #13) with differing sign of Am^. It should 
be noticed that this is true no matter how large the matter effect, quite contrary to the 
case of CP measurement. Therefore, determination of the signs of Am^ is impossible in a 
single T measurement experiment unless one of the following conditions is met; (a) one of 
the solutions is excluded, for example, by the CHOOZ constraint, or (b) some additional 
information, such as energy distribution of the appearance electrons, is added. 

We emphasize that a complete understanding of the structure of the parameter degen- 
eracy should be helpful for one who want to pursue solution of the ambiguity problem in an 
experimentally realistic setting. We, however, do not attempt to discuss the #23 <-> f — #23 
degeneracy. We also do not try to solve the problem of parameter degeneracy exactly though 
it is in principle possible by using an exact but reasonably compact expression of the oscilla- 
tion probability obtained by Kimura, Takamura, and Yokomakura ||20|| . Instead, we restrict 
ourselves into the treatment with the approximation introduced by Burguet-Castell et al. 
H in which the approximate formula for the oscillation probability derived by Cervera et 
al. |21] was employed. Though not exact, it gives us much more transparent overview of the 
problem of parameter degeneracy. 
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II. PROBLEM OF PARAMETER AMBIGUITY IN CP AND T VIOLATION 

MEASUREMENT 



We define the "CP (T) parameter ambiguity" as the problem of having multiple solutions 
of (5, 9i 3 ) and the sign of Am^, for a given set of measured values of oscillation probabilities 
of P(i/„ -> v e ) and its CP (T) conjugate, GP[P(u)] = P(P A1 -> z> e ) (T[P(z/)] = P(z/ e -> ^)). 
We concentrate in this paper on this channel because precise measurement is much harder 
in other channels, e.g., in v e — * z/ r . Our use of — > u e and its CP-conjugate is due to 
our primary concern on conventional super-beam type experiments [22|. The reader should 



keep this difference in mind if they try to compare our equations with those in Refs. P|T2" 



in which they use v e — > and its CP-conjugate, a natural choice for neutrino factories 
pl|,|23| . It should also be noted that the neutrino factories and the superbeam experiments 
are studying processes which are T-conjugates. 



In this section we utilize the CP and the T trajectory diagrams introduced in JT0[ and 
||12|| , respectively, to explain what is the problem of parameter ambiguity and to achieve 
qualitative understanding of the solutions without using equations. But before entering into 
the business we want to justify, at least partly, our setting, i.e., prior determination of all 
the remaining parameters besides 5 and #13. 

A. Problem of parameter degeneracy; set up of the problem 

We assume that at the time that an experiment for measuring (5, 813) is carried out all 
the remaining parameters in the MNS matrix, 823, |Am| 3 |, 812, and Am^, are determined 
accurately. The discussion on how the experimental uncertainties affect the problem of 
parameter degeneracy is important, but is beyond the scope of this paper. It should be 
more or less true, because 823 and Am| 3 will be determined quite accurately by the future 



long-baseline experiments [ 24 - 26 1 up to the 823 <-> f — 6*23 ambiguity. Most notably, the 



expected sensitivities in the JHF-SK experiment in its phase I is [24 
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5(sin 2 2 9 23 ) ~ 10~ 2 

5(|Am 2 3 |) ~ 4 x 10~ 4 eV 2 (1) 

at around |Am 23 | = 3 x 10~ 3 eV 2 . On the other hand, the best place for accurate determi- 
nation of 6\2 and Am 2 2 , assuming the LMA MSW solar neutrino solution, is most probably 
the KamLAND experiment; the expected sensitivities are f2"7fl 

5(tan 2 12 ) ~ 10% 

<5(Am 2 2 ) ~ 10% (2) 

at around the LMA best fit parameters. Therefore, we feel that our setting, prior determi- 
nation of all the mixing angles and Am 2, s besides 5, #i 3 and the sign of Am 23 , is reasonable 
one at least in the first approximation. 



B. Pictorial representation of parameter ambiguities 



In this subsection we use CP and the T trajectory diagrams |jnjl2] to explain intuitively 
what is the problem of parameter degeneracy, and to achieve qualitative understanding of 
the solutions without using equations. In Fig. 1 we display four CP trajectories in the 
P-CP[P] bi-probability space which all have intersection at P(v) = 1.9% and CP[P(V)] = 
2.6%. The four CP trajectories are drawn with four different values of # 13 , sin 2 2# 13 = 0.055, 
0.050, 0.586, 0.472, and the former (latter) two trajectories correspond to positive (negative) 
Am 2 3 . The analytic expressions of the four degenerate solutions will be derived in Sec. IV. 
The neutrino energy E and the baseline length L are taken as E = 250 MeV and L = 130 
km, respectively^]. The setting anticipates an application to the CERN-Frejus project [28 



where the regions with parameter ambiguities are widest. The values of all the remaining 
mixing parameters are given in the caption of Fig. 1. 



In all the figures of this paper we do not average over neutrino energy distributions as was 
performed in Refs. [10,12| but use a fixed neutrino energy specified for each figure. 
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FIGURES 
L = 130 km, E = 250 MeV 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

P(v) (%) 

FIG. 1. An example of the degenerate solutions for the CERN-Frejius project in the 
P{y) = P{V)jl — * Ve) verses CP[P{v)] = P{pu — ¥ v e ) plane. Between the solid (dashed) lines 
is the allowed region for positive (negative) Am 2 3 and the shaded region is where solution for 
both signs are allowed. The solid (dashed) ellipses are for positive (negative) Am 2 3 and they all 
meet at a single point. This is the CP parameter degeneracy problem. We have used a fixed 
neutrino energy of 250 MeV and a baseline of 130 km. The mixing parameters are fixed to be 
|Am 2 3 | = 3 x 10- 3 eF 2 , sin 2 2<9 23 = 1.0, Amf 2 = +5 x 10" 5 eF 2 , sin 2 29 l2 = 0.8 and Y e p = 1.5 g 
cm -3 . 
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Figure 1 demonstrates that there is a four-hold degeneracy in the determination of (5, 
#13) for a given set of P(v) and CP[P(z/)]. The region between the solid lines and the region 
between the dashed lines in Fig. 1 are the regions where two-fold degeneracies exist in the 
solutions of (5, #13) for positive and negative Am^ sectors, respectively. (See Eq. ( pp in 
Sec. IV.) It is intuitively understandable that the region where degenerate solutions exists is 
the region swept over by the CP trajectories when the parameter # 13 is varied, while keeping 
other mixing parameters and the experimental conditions fixed. The shaded region is the 
region where the full four-hold degeneracy exists. 

Now we turn to the T measurement. In Fig. 2 we display four T trajectories in the P- 
T[P] bi-probability space which all have intersection at P(v) = 1.7% and T[P(z/)] = 2.5%. 
The four T trajectories are drawn with four different values of #13, sin 2 2# 13 = 0.05, 0.0427, 
0.575, 0.490, and the former (latter) two trajectories correspond to positive (negative) Am^. 
Matter effects split the positive and negative AmJ 3 trajectories, see [12], thus different values 
of 9 are required for mixed sign trajectories to overlap. The remaining mixing parameters 
and the experimental setting are the same as in Fig. 1. 
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L = 130 km, E = 250 MeV 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

P(v) (%) 



FIG. 2. An example of the degenerate solutions using the energy and path length of the 
CERN-Frejius project in the P{y) = P(vn —* v e ) verses T[P(v)] = P{y e — ► Vn) plane. For this 
figure there is complete overlap in the region (shaded) that allows solutions for either sign of Amf 3 . 
The solid (dashed) ellipses are for positive (negative) Amf 3 and they all meet at the "measured 
point", (P, P T ) = (1.7, 2.5)%. This is the T-parameter degeneracy problem. Notice that for the 
each ellipse with positive Amj 3 there is a completely degenerate ellipse with negative Amf 3 . This 
feature will be explained in Sec. IIIc. Parameters are the same as in Fig. El 
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A clear and interesting difference from the CP diagram manifests itself already at this 
level, reflecting the highly symmetric nature of the T conjugate probabilities as will be made 
explicit in eq. @. Namely, the two different- Am^-sign diagrams (the first and the third, and 
the second and fourth) completely overlap with each other. In the next section we will make 
it clear that the complete degeneracy originates from a symmetry. Therefore, discrimination 
of the sign of Am^ is impossible in a single T-violation measurement experiment unless one 
of the solutions is excluded by an other experiment. We will demonstrate in Sec. Ill that the 
degeneracy is not accidental one specific to this particular case, but its existence is generic. 
There is always a four-hold degeneracy in T measurement. 



III. PARAMETER DEGENERACY IN T-VIOLATION MEASUREMENTS 

We start by presenting an analytic treatment of the problem of parameter degeneracy 
in T-violation measurements primarily because it is simpler and instructive. To do this 
we generalize the formalism developed by Burguet-Castell et al. || by treating the cases 
of positive and negative Am 2 3 simultaneously. It provides basis of our treatment of the 
degeneracy between the same as well as across the alternating Am^-sign probabilities. 
It will become clear from the following discussions that the treatment of the mixed-sign 
degenerate solutions, for both CP and T measurements, can be done as a straightforward 
generalization of the same-sign degeneracy case by simply taking account of duplication due 
to the alternating sign of Am^ [10 . 



There are four basic equations satisfied by the T-conjugate probabilities in the case of T 
measurement for small sin #13 : 

p(u)+ = x + e 2 + Y + e cos (5 + + p 

T[P(u)) + = P» + = X + 6 2 + Y + 6 cos (5 - + P & 

P( v )_ = X„6 2 + Y_6 cos (5 - + Pq 

T[P{y)\- = P T (v)- = X_6 2 + Y_6 cos (s + ^) + P & (3) 
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where X± and Y± are given in Appendix, P G indicates the term which is related with solar 



neutrino oscillations, and A 



lAmfJL 



13 



2E 



Note that ± here refers to the sign of Am 2 3 and 8 



is an abbreviation of #13 ~ sin #13. In the next subsections we discuss the possible solutions 
for 8 and S for a given measurement of both P and P T for both positive and negative sign 
of Am 2 3 . 

A. The same-sign degeneracy; T measurement 

The treatment in this subsection applies for two overlapping T trajectories with the same 
sign of Amf 3 . The degeneracy associated with alternating-sign trajectories will be explored 
in the next subsection. 

There are two sets of approximate solutions of (£J), 9% and S\ , where (i = 1,2) and (i = 
3,4) denotes the solutions in the positive and negative Am^ sectors, respectively. They are 

A 13 ' 



8; 



IP-P<* 



Y 



± 



X. 



± 



2X4 



COS Si ± 



P T -Pc, 



Y. 



± 



cos Si T 



2 
A 13 



(4) 



1 \\ X± 2X ± ' ~ V 1 1 2 
where ± correspond to solutions in positive and negative Amf 3 sectors 0. We then obtain, 



e.g., for the positive Am^ sector 



8 2 — 81 
82 — 8\ 



2X + 

y+ 

2X, 



J^COS^ — COsJi) COS 
(^cos52 — COS(5i) COS 



A 



13 



A 



13 



— (sin S2 — sin Si) sin 
+ (sin S2 — sin Si) sin 



A 



13 



2 

A13 



(5) 



which entails the degeneracy that if (81, Si) is a solution so is 

82 = 81 + — — cos Si cos f — — ^ and 5 2 = it — Si, 
X 4. V 2 / 



(6) 



2 The above solutions are exact solutions to the system of Eq. (^) if we were to add terms 
Y± cos 2 ± % a ) /(4X±) to the equations, @ . In what follows we have systematically ignored 
terms of 0(Y£/X±). 
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in addition to the trivial solution. A similar degeneracy holds also for the negative Amf 3 
sector, that is, if (#3, 83) is a solution so is 

#4 = #3 + cos 83 cos an d 84 = n — 83. (7) 

Both of these same sign Atoj 3 degeneracies are in matter though they look like the vacuum 
degeneracies as discussed in |[10|| . Notice that if the experimental setup is chosen such that 
cos f ^f 3 -) = ||] or nature has chosen cos 8 = then the same sign degeneracies are removed. 



B. The mixed-sign degeneracy; T measurement 



Let us now examine the problem of parameter degeneracy which involves positive and 
negative Amf 3 . The basic equations fl3|) can be approximately solved for mixed sign situation 
as: 

A 13 ' 



0i 
9, 



IP-R 







X 4 



2X+ 



'P-Pe 



X_ 



P T -P, 







2X 
2X7 



cos \ 5i + 

- COS I 8\ 



cos ( 8 



cos ( 83 + 



2 

A13 

2 

¥) 

A 13 



(8) 



\ X. 2X_ ' " V' 2 
We will now exactly determine 9 and 5 using the above set of approximate solutions 
and (U), as our starting point. First, 

Ai 3 VX7 



(9) 



sin 5 1 sin 
sin £3 sin 



2 

A13 



Y + 
fXl 

~y7 



P-Pc 







p t -p 







P-Po, 



P T -P P 



(10) 

(11) 



and cos5j is given by cos5j = ±y 1 — sin <5j. Using these cos5j the values of 9 are given by 

A13 



9i 



VP -P Q + JP T -P £ 



2^fX~Z 
2~JX = 



Y+ 
2X H 

r_ 

2X1 



COS 5i COS 



COS 83 cos 



2 

A13 



(12) 
(13) 
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To relate these alternating sign solutions we use the identity, see Appendix, 

xfXl xfXZ 



(14) 



Y+ Y- 

derivable under the Cervera et al. approximation^. Then, it follows that 

sin 5i = sin S 3 

(cos5i + cos5 3 ) cos ^ = -2^± (Jx^9 1 - JxZe 3 ) (15) 

One can choose without loss of generality S 3 = n — Si as a solution of (|i~5l) . Then, for a given 
P and P T measurement, apart from (9i, Si) there are three other solutions^ given by 

9 2 = 9i + cos 5\ cos f — r— and 5 2 = n — 5i 



^=\j§ te i and S 3 = ir-5 1 (16) 

4 = 9 3 — -^r- cos Si cos ^— and 5 4 = 5i 

Therefore, there is no ambiguity in determination of 5 in T violation measurement apart 
from the one 5 — ► 7r — 5 independent of the sign of Amf 3 . Fortunately, this degeneracy does 
not obscure existence or non-existence of leptonic T (or CP) violation. This feature arises 
because of highly constrained nature of system (|3|) of T-conjugate probabilities. 

The physically allowed region of the T diagram is determined by the constraint that 
sin 2 Si < 1 which in terms of P and P T is 



X , V 2 / X_ V 2 



3 Unless Eq. (|14|) holds we get into trouble because then Eqs. (||) or (|9|) does not allow the (same- 



sign) solution 9\ = 62 and Si = 82, which must exist as shown in Ref. [12]. Therefore, use of the 
formula of oscillation probability obtained by Cervera et al. who summed up all order matter effect 
is essential. 

As a convention, we have chosen cos 5i < so that for A13 < tt, 61 > 62 and 9 3 > 64. 
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and is the same region for both signs of Am^ because of the identity Eq. (|T4"D . In Fig. |2|, 
this region is the shaded region using the CERN-Frejius parameters. At the boundary of 
the allowed physical region cos 5 = and the same sign degeneracy vanishes, however, the 
opposite sign degeneracy is non-zero. 

We define the fractional differences, A9/9, by 



A9 



Oj — On 



(18) 



to quantify how different the two degenerate solutions are. In fact, one can obtain simple 
expressions for the various fractional differences; 

-2Y + cos S x cos (^f ) 



A9 



A9 



'A£ 



A£ 



12 



31 



11 



32 



'A£ 



34 



12 



Oci(yp=p^ 
UfXl - JXD 



P T -P & 



on 



or 



-2Y + cos^i cos (4 



OC 



or 



^[yr^rPo + y/FT -p & 

—2Y + cos S\ cos 



2 



P T -P & 



+ 2 



o^ + yxz 

OCZ-^fxO 



on 



or 



(19) 

(20) 
(21) 

(22) 



The same sign fractional difference, (1,2) and (3,4), decreases with increasing P and P T and 
thus 9, whereas the first mixed sign fractional difference, (3,1) and (4,2), is independent of 
the size of P and P T and thus 9. The second and third mixed sign fractional differences, 
(1,4) and (2,3), are similar to the same sign fractional difference but off set by an energy 
dependent constant. The relationship between the fractional difference in the measured 
quantity sin 2 29 and 9 is simply 



sin 2 20 i - sin 2 29 ' ■ 



9j — 9 a 
2 3 



+ 9^/2 =2 \T / 



(23) 



(sin 2 29 i + sin 2 2%)/2 

for 1 > 9i,9j > \9,i -9j\. 

In Fig. 3 thru 5 we have plotted the differences in the allowed 9 solutions divided by 
half the sum for the CERN-Frejus, JHF-SK and FNAL-NuMI [^9j possible experiments 
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using — > z/ e and its T-conjugate f e — > v^. The regions where this fractional difference 
is small are regions where the parameter degeneracy inherent in such measurements is only 
important once the experimental resolution on 9 for a fixed solution is of the same size or 
smaller. Notice that near the boundaries on the allowed region the fractional differences are 
small for the same sign solutions. For the mixed sign, either the (1,4) or (3,2) fractional 
difference plots have a line for which the fractional difference is zero. This line can be 
understood as follows; for a given small value of 9 the positive and negative Am^ ellipses 
overlap and intersect at two points. As 9 varies these intersection points give us this line 
with zero mixed sign fractional difference. 



A0/9 (%) for L =130 km, E = 250 MeV 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

P(V) (%) 

FIG. 3. The iso-fractional differences, as a %, for the allowed solutions for the mixing angle #13 
in the P{v) = P{Vfj, — ► Ve) verses T\P{v)\ = P{v e — ► v^) plane for the CERN-Frejius project. The 
fractional differences for solutions (3,4) is identical to that for (1,2) and the fractional difference 
for (3,2) equals the fractional difference for (1,4) plus or minus a constant, see eq. ©-(H). In this 
case the fractional difference (3,2) has a zero contour. The parameters are the same as in Fig. |]. 
The ellipses are labelled to show the relevant size of sin 2 2#i3 for this figure. 
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A0/9 (%) for L =295 km, E = 1 GeV 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

P(V) (%) 

FIG. 4. The iso- fractional differences, as a %, for the allowed solutions for the mixing angle 
#13 in the P{y) = -P(^7t - ► Ve) verses T[P{v)\ = P{y e — > v^) plane for the JHK-SK project. The 
fractional differences for solutions (3,4) is identical to that for (1,2) and the fractional difference 
for (3,2) equals the fractional difference for (1,4) plus or minus a constant, see eq. (p~9|)-(p2|) . The 
zero contour appearing in the (1,4) fractional difference is explained in the text. The parameters 
are the same as in Fig. ||. The ellipses are labelled as in Fig.||. 
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A0/9 (%) for L =732 km, E = 2 GeV 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

P(V) (%) 

FIG. 5. The iso-fractional differences, as a %, for the allowed solutions for the mixing angle #13 
in the P(y) = P(> M -> v e ) verses T[P(v)} = P(v e -> v^) plane for the FNAL-NUMI project. The 
fractional differences for solutions (3,4) is identical to that for (1,2) and the fractional difference 
for (3,2) equals the fractional difference for (1,4) plus or minus a constant, see eq. (|T9|)-(|22|), At 
very small probability the (1,4) fractional difference has a zero contour. The parameters are the 
same as in Fig. H. The ellipses are labelled as in Fig.||. 

C. Symmetry between the two alternating-Amf 3 -sign T diagrams 

The observant reader will notice that there is in a general one-to-one correspondence 
between the solutions with positive Am^, labelled 1 and 2, and those solutions with negative 
Am? 3 , labelled 3 and 4, in Eq. ([16]) by using the identity, Eq. (|T4]). In fact this correspondence 
applies not only to the solutions but to the complete T diagram. For a given T trajectory 
with positive Am^ there always exists a T trajectory with negative Am^ with a different 
value of 8, which nevertheless completely overlaps the positive trajectory, see Fig. [| This 
surprising phenomenon occurs because there exists a symmetry in T probability system 
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defined in Eq. (|3]). 

One can observe from eq. (|3|) that a positive Am^ trajectory which is defined by the first 
two equations of @ can be used to generate a negative Am^ trajectory which completely 
overlaps with the original one by the transformation 

5 —> 7T — 5 



,~fco (24) 

This means that for a measure set of P and P T there is always two set of solutions with 
different sign of Am^. There is no way to resolve this ambiguity because the two T trajec- 
tories are completely degenerate. It should be noticed that this situation occurs no matter 
how large a matter effect at much longer baseline. In such a case, two T trajectories with the 
same #13 but opposite sign of AmJ 3 are far apart, and one would have expected that there 
is no ambiguity in determination of the sign of Amf 3 . Hence, there is a "no-go theorem" for 
the determination of the sign of Amj 3 by a single T violation measurement. The possible 
cases in which the "theorem" is circumvented are, as mentioned in Introduction, (a) one of 
the solutions is excluded, for example, by the CHOOZ constraint, or (b) some additional 
information, such as energy distribution of the appearance electrons or another T-violation 
measurement with different parameters, is added. 

IV. PARAMETER DEGENERACY IN CP- VIOLATION MEASUREMENTS 

We now turn to the analytic treatment of the parameter degeneracy in CP-violation 
measurements. We proceed in an analogous way to the analytic treatment of T-violation 
given in sec. III. 

We start with the four basic CP equations for small sin# 13 : 

p(u)+ = x + e 2 + Y + e cos (5 + + p q 
cp[p(v)} + = P(u)+ = x + e 2 + F + flcos (s-^j + P 
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P(v)_ = X_6 2 + Y_9 cos ( 5 



A 



13 



CP[P(i/)]_ =P(^)- = X_# 2 + F„flcos (<J 



+ P Q 

Aia 



P 







(25) 



where X± and Y± are given in Appendix. As before, P Q indicates the term which is related 
with solar neutrino oscillations, A 13 = — the ± here refers to the sign of AmJ 3 and 
is an abbreviation of 9 13 ~ S13. 

Note that there exist relations among coefficients; 



X± = X(±Am 2 3 ,a) 
X± = X(±Am 2 3 , -a) 



(26) 
(27) 



In leading order in , there exist further relations, 



X , = X_ and X_ = X, 



(28) 



which follows from the CP-CP relation [|T^] (see Appendix). Finally, it follows under the 
approximation of Cervera et al. [2~I| that 



y, = -y_, r_ 



(29) 



We fully utilize the symmetry relationships ( |28| ) and (|29|) in the unified treatment of 
the same-sign and the mixed-sign degeneracies. The basic equations ( p5[) can be solved for 
generic mixed sign situation as: 



0i 
0i 

03 
03 



/P - P^ F, 



X 4 



2X, 



cos ( Si 



A 



13 



P-Pm k 



+ 2x: cos(51 



IP - p m 



X_ 



2X_ 



cos ( 5; 



P-Pn Y 



+ 



^ X+ ^2X 4 
The solution of these equations are: 



cos 5, 



2 

A13 

2 

A13 

2 
A 13 



(30) 
(31) 
(32) 
(33) 
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1*1,2= n 



1 


-C {+) sin 


A13 


D 




2 


1 


-C (+) sin 


A13 


D 




2 



Ai 



#1,2 — 



1 

2DH 



\P - P )X 4 



~ Pe)X- . x . A 13 
h siik)i )2 sin 



#3,4 -^py 



-P & )X_ ^(P - P Q )X + . . Ax 
-yT + sm J 3,4 sin — 



13 



(34) 
(35) 

(36) 
(37) 



where 



D^ C M 2 sin 2 (^) + CH 2 cos 2 (^p) 



2 ^ X_ 



2 \ n 



A 

Y. 
X 

X_ 



AP + = 



/P-P 



X 4 



P-P, 



(38) 
(39) 

(40) 
(41) 



The sign in (j35|) is determined relative to ( 34 ) so that it reproduces the pair of degenerate 
solutions in the case of a precisely determined value of 9 ri |TJJ. It should be noticed that 



provided \j D — (AP±) 2 is real the constraint | sin ^ | < 1 is satisfied automatically in (|34"D 
and (H). 

Let us focus first on the features of the same-sign degenerate solution. The set (0 i: 5i) 
with i = 1,2 (z = 3, 4) describes two degenerate solutions with positive (negative) Am^ for 
given values of P and P. Of course, they reproduce the relationships obtained by Burguet- 
Castell et al in Pi : 



sin <5 2 — sin Si = —2 



9 2 — 0i = 



'sin #i + z cos 5i 
, 1 + z 2 



sin Sx + z cos 5i \ 

C(+) 2 _ 



sin 



^13 



(42) 
(43) 



where 



CH tan 



(44) 
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Let us illuminate how the relative phases between <5's between these degenerate solutions 
can be obtained in a transparent way. Toward the goal we first calculate cos 5;. cos 5^2 and 
cos $3,4 can be obtained from (|34|) and (35) , respectively, by replacing of by =fC ( - =f - ) and 



sin 4p by cos ^ a and vice versa. One can show by using these results that 



2 J 2 

1-z 2 
l + z 

which implies that 



cos (Si + 8 2 ) = cos (S 3 + 5 4 ) = n , _ 2 , (45) 



S 2 = 7T - Si + arccos((^ 2 - l)/(z 2 + 1)) 

5 4 = 7T - 5 3 + arccos((^ 2 - l)/(z 2 + 1)). (46) 

Thus in the allowed region of bi-probability space S 2 (64) differs from it — Si (tt — S3) by 
a constant, arccos((z 2 — l)/(<2 2 + 1)), which depends on the energy and path length of the 
neutrino beam but not on the mixing angle 9. Near the oscillation maximum, z — > 00, this 
constant vanishes so that S 2 — 7r — Si and £3 ~ 71 — S4 as noticed in |5| . 
For the mixed-sign degenerate solution one can show that 



( AP + AP_ - JD - APl^/D - AP. 



D 



cos (Si - S 3 ) = 

■ , x x , (AP +v /P-AP 2 + AP-y/P-AP 2 
sm (S 1 - S 3 ) = ± i s (47) 

One can show, for example, cos (Si — S3) = —1 and sin (Si — S3) = in the P — > P limit by 
noting that AP_ = — AP + in the limit. It means that £3 = Si + 7r (mod. 27r), in agreement 



with the result obtained in Ref. [12 



The conditions for existence of the same-sign solution are 

D - (AP+) 2 > and D- (AP_) 2 > (48) 

for positive and negative Am 2 3 , respectively. The condition for existence of the mixed- 
sign solution is the intersection of the two regions which satisfy the conditions of eq. (^8|) . 
An example of the regions satisfying conditions for existence of the same-sign as well as 
mixed-sign solutions are depicted in Fig. 1. 
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The maximum value of P and P which allows mixed sign solutions is determined by 
D — (AP + ) 2 = D — (AP_) 2 = with P = P. This occurs for a critical value of P given by 

X + X_D , . 

p - = p ° + wv^y (49) 

which can be used to determine the critical value of 9 as 

C(+) s in 2 ^ , 

s ~ = ^7W- (50) 

There is no degeneracy in the value of 9 at this critical point, i.e. 9\ = 92 = #3 = 0&- 
An example of this can be seen in Fig. || At the first peak in the oscillation probability, 
A 13 = 7T, the value of the critical 9 is simply given by 

WA» = »)= !%!-— ■ (51) 



As A13 — > 27T, the critical goes to zero and so does the oscillation probabilities P and P 
as this is the position of the first trough in the oscillation probabilities. 

In Fig. 6 thru 8 we have plotted the fractional difference, (^r), see eq. ([18]), for the 
CERN-Frejius, JHF-SK and FNAL-NUMI possible experiments using — > u e and its CP- 
conjugate ^ v e . The regions where this fractional difference is small are regions where 
the parameter degeneracy inherent in such measurements is only important once the ex- 
perimental resolution on 9 for a fixed solution is of the same size or smaller. Notice that 
near the boundaries on the allowed region the fractional differences are small for the same 
sign solutions. For the mixed sign, (1,3), fractional difference plots there is a zero along the 
diagonal. This is explained by the fact that in our approximation the positive and negative 
Am 2 3 ellipses, for a given 9, intersect along the diagonal. Fig. 9 thru 11 are similar to 
Fig. 6 thru 8 expect that the size of solar Am 2 2 has been increased by a factor of two to 
1 x 10 _4 eV 2 . Notice that the parameter degeneracies problem is more pronounced as the 
solar Am 2 is increased. At even larger values of Am 2 2 , our approximations become less 
reliable. 



23 



A0/9 (%) for L =130 km, E = 250 MeV 
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FIG. 6. The iso-fractional differences, as a %, for the allowed solutions for the mixing angle 
#13 in the P{y) = P {ya —* v e) verses CP\P(y)\ = P{p^ — > v?) plane for the CERN-Frejius project. 
The parameters are the same as in Fig. [l] with Am 2 2 = 5 x 10~ 5 eV 2 . The dashed ellipses are 

labelled CP'.^'. to show the relevant size of sin 2 for this figure. 

sir 26»i3 
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A9/9 (%) for L =295 km, E = 1 GeV 
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FIG. 7. The iso- fractional differences, as a %, for the allowed solutions for the mixing angle 
#13 in the P{y) = -P(^ - ► v e ) verses CP[P(y)\ = P{p^ — »■ ^ e ) plane for the JHF-SK project. The 
parameters are the same as in Fig. [l] with Amf 2 = 5 x 10 -5 eV 2 . The ellipses are labelled as in 



Fig§ 
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A0/0 (%) for L =732 km, E = 2 GeV 
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FIG. 8. The iso-fractional differences, as a %, for the allowed solutions for the mixing angle #13 
in the P{y) = Piy^ u e ) verses CP[P(v)} = P{v^ u e ) plane for the FNAL-NUMI project. The 
mixed sign fractional differences, (A9/8) l3 and 14 terminate at around P = CP[P] « 1.6% because 
above this probability the sign of Am^ is determined, as discussed in the text. The critical value 
of P (~ 1.6%) and sin 2 26\z (~ 0.033) can be calculated from eq. (f49|) and (^0|). The parameters 



are the same as in Fig. [l] with Am 2 2 = 5 x 10 5 eV 2 . The ellipses are labelled as in Fig|| 
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A0/0 (%) for L =130 km, E = 250 MeV 
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FIG. 9. The iso-fractional differences, as a %, for the allowed solutions for the mixing angle 
#13 in the P{y) = P(yu —* v e) verses CP\P(y)\ = P{p^ — > P e ) plane for the CERN-Frejius project. 
The parameters are the same as in Fig. [I] except for the solar Amf 2 which is set to 1 x W~ eV 2 
for this plot. The ellipses are labelled as in Fig|| 
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A0/0 (%) for L =295 km, E = 1 GeV 
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FIG. 10. The iso-fractional differences, as a %, for the allowed solutions for the mixing angle 
#13 in the P{y) = P{yu — ► v e ) verses CP[P(y)\ = P{p^ — »■ ^ e ) plane for the JHF-SK project. The 
parameters are the same as in Fig. [I] except for the solar Amf 2 which is set to 1 x 10~ 4 eV 2 for this 
plot. The ellipses are labelled as in Fig|| 
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A0/0 (%) for L =732 km, E = 2 GeV 
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FIG. 11. The iso-fractional differences, as a %, for the allowed solutions for the mixing angle 
13 in the P(v) = Pfa f e ) verses CP[P(v)] = Pfa -> P e ) plane for the FNAL-NUMI project. 
The parameters are the same as in Fig. [I] except for the solar Amf 2 which is set to 1 x W~ eV 2 
for this plot. The ellipses are labelled as in Fig|| 

V. SUMMARY AND CONCLUSION 

In this paper we have given a complete analytic treatment of the parameter degeneracy 
issue for #13, sign of Am^ and the CP and T violating phase 5 that appears in neutrino 
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oscillations. For a given neutrino flavor transition probability and its CP or T conjugate 
probability we have derived the allowed values of #13, sign of AmJ 3 and the CP and T 
violating phase 5. We have given explicit expressions of degenerate solutions and obtained, 
among other things, exact formulas for the relationship between solutions of 5's up to the 



In general there is a four-fold degeneracy, two allowed values of #13 for both signs of 
Aml 3 . This is always true for the T violation measurement whereas for the CP violation 
measurement the four-fold degeneracy can be reduced to two-fold degeneracy if matter ef- 
fects are sufficiently large, or we live close to the region 5 ~ 7r/2 or 3n/2. The significance of 
matter effects dependence on the energy of the neutrino beam, the separation between the 
source and the detector as well as the density of matter between them. The fractional differ- 
ence of 6*i3 between the various solutions has been calculated which can be compared with 
the experimental sensitivity for a given setup to determine whether or not the degeneracy 
issue is significant or not. 

For the possible future experimental setups CERN-Frejius, JHK-SK and FNAL-NUMI 
we have given numerical results for the channel — > v e and its CP and T conjugate. The 
CP conjugate being most relevant for these future Super-beam experiments. For the CERN- 
Frejius, JHF-SK and FNAL-NUMI experimental setups the parameter degeneracy issue is 
only relevant once the experimental resolution on the determination of #13 is better than 
15%, 10% and 5% respectively, assuming a transition probability near 1% and a Am^ 2 = 
5 x 10~ 5 eV 2 , see Fig. || - |8| . At larger values of Am^ the parameter degeneracy issue 
becomes more important. These iso-fractional difference plots are useful for comparing the 
sensitivity of different experimental setups, neutrino energy, path length and experimental 
sensitivity, to this parameter degeneracy issue. 
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VI. APPENDIX 

The standard flavor transition probability for neutrino oscillations in the — > v e channel 
can be written as 

P{y)± = X±9 2 + Y ± 9 cos (s ± -^j + P Q (52) 

where the ± signs in X± and Y± refer to positive or negative values of Am^, 9 is an 
abbreviation for sin # 13 and P indicates the terms related to solar neutrino oscillations. For 
details on the approximations used in deriving this transition probability see ref. pi| . All 



other channels used in this paper, u e — * and —>■ u e , can also be expressed with the 
same variables, see Sec. Ill and IV. 

The coefficients X± and Y± are determined by 
,2 ( ^ri\ . 2 ( B t 



X ± = 4sl 3 ^-j sin^^tj, (53) 



Y ± = ±8c 12Sl2 c 23S2 3 ^ ) ) sin i^—j sin {—) (54) 

P = 4sin 2 2^ 12 (^) 2 sin 2 (^) (55) 

with 

\AmfAL , , 

Ay = o ; j ' and B ± = \aL ± A 13 |, (56) 

where a = \f2GpN e denotes the index of refraction in matter with Gp being the Fermi 
constant and N e a constant electron number density in the earth. 
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Obviously from the above definitions, X± and Y± satisfy the identity 




which is used throughout this paper. 



32 



REFERENCES 



[1] Kamiokande Collaboration, Y. Fukuda et al, Phys. Lett. B335 (1994) 237; 

Super-Kamiokande Collaboration, Y. Fukuda et al, Phys. Rev. Lett. 81 (1998) 1562; 
ibid. 85 (2000) 3999. 

[2] Homestake Collaboration, K. Lande et al, Astrophys. J. 496 (1998) 505; 

SAGE Collaboration, J. N. Abdurashitov et al, Phys. Rev. C 60 (1999) 055801; 
GALLEX Collaboration, W. Hampel et al, Phys. Lett. B 447 (1999) 127; 
Super-Kamiokande Collaboration, S. Fukuda et al, Phys. Rev. Lett. 86 (2001) 5651; 
ibid. 86 (2001) 5656; 

SNO Collaboration, Q. R. Ahmad et al, Phys. Rev. Lett. 87 (2001) 071301; ibid. 89 
(2002) 011301; 89 (2002) 011302. 

[3] K2K Collaboration, S. H. Ahn et al, Phys. Lett. B 511 (2001) 178; 

K. Nishikawa, Talk presented at XXth International Conference on Neutrino Physics 
and Astrophysics (Neutrino 2002), May 25-30, 2002, Munich, Germany. 

[4] Z. Maki, M. Nakagawa and S. Sakata, Prog. Theor. Phys. 28 (1962) 870. 



[5] M. B. Smy, |hep-ex/020202D| ; SNO Collaboration, Q. R. Ahmad et al, 89 (2002) 011302; 



J. N. Bahcall, M. C. Gonzalez-Garcia, and C. Pena-Garay, |hep-ph/0204314| ; V. Barger 
et al., Phys. Lett. B 537 (2002) 179; A. Bandyopadhyay et al, |hep-ph/0204286| ; 



PC. de Holanda, AYu. Smirnov |liep-ph/0205"24ll ; G. L. Fogli et al, |hep-ph/ 0206162 



M. Maltoni et al, |hep-ph/ 020 7227 . 



[6] S. P. Mikheev and A. Y. Smirnov, Nuovo Cim. C 9 (1986) 17; 
L. Wolfenstein, Phys. Rev. D 17 (1978) 2369. 

[7] CHOOZ Collaboration, M. Apollonio et al, Phys. Lett. B 420 (1998) 397; 

ibid. B 466 (1999) 415; See also, The Palo Verde Collaboration, F. Boehm et al, Phys. 
Rev. D 62 (2000) 072002. 

33 



[8] T. Kajita, H. Minakata, and H. Nunokawa, Phys. Lett. B528 (2002) 245 flnep- 
ph/01 12345| . 



[9] J. Burguet-Castell, M.B. Gavela, J.J. Gomez-Cadenas, P. Hernandez and O. Mena, 
Nucl. Phys. B 608 (2001) 301. 



[10] H. Minakata and H. Nunokawa, JHEP 0110 (2001) 001 ||hep-ph/ 01 08085)1 . 



[11] V. Barger, D. Marfatia, and K. Whisnant, Phys. Rev. D 65 (2002) 073023. 

[12] H. Minakata, H. Nunokawa, and S. J. Parke, Phys. Lett B, B537 (2002) 249 Eepj 



ph/0204171|| . 



[13] J. Arafune and J. Sato, Phys. Rev. D55 (1997) 1653; J. Arafune, M. Koike and J. Sato, 
Phys. Rev. D56 (1997) 3093 [Erratum ibid. D 60 (1999) 119905]; 
H. Minakata and H. Nunokawa, Phys. Rev. D57 (1998) 4403; Phys. Lett. B413 (1997) 
369; ibid. B495 (2000) 369; 

K. Dick, M. Freund, M. Lindner and A. Romanino, Nucl. Phys. B 562 (1999) 29; 

O. Yasuda, Acta. Phys. Polon. B 30 (1999) 3089; M. Koike and J. Sato, Phys. Rev. 

D61 (2000) 073012; Erratum ibid. D 62 (2000) 079903. 

[14] S. J. Parke and T. J. Weiler, Phys. Lett. B501 (2001) 106. 



[15] H. Minakata and H. Nunokawa, Nucl. Phys. Proc. Suppl. 110 (2002) 404 flhep 



ph/0111131 



[16] V. Barger, D. Marfatia, and K. Whisnant, |hep-ph/ 0206038 



[17] P. Huber, M. Lindner, and W. Winter, |hep-ph/0204352 . 



[18] A. Donini, D. Meloni, and P. Migliozzi, |hep-ph/0206034 



[19] J. Burguet-Castell, M.B. Gavela, J.J. Gomez-Cadenas, P. Hernandez and O. Mena, 
|hep-ph/020708q . 



34 



[20] K. Kimura, A. Takamura and H. Yokomakura, Phys. Lett. B537 (2002) 86; |licp 



ph/0205295 . 



[21] A. Cervera, A. Donini, M. B. Gavela, J. J. Gomez Cadenas, P. Hernandez, O. Mena 
and S. Rigolin, Nucl. Phys. B579 (2000) 17 [Erratum-ibid. B593 (2000) 731]. 

[22] The idea of low-energy conventional "superbeam" for measuring leptonic CP violation 
may be traced back to 

H. Minakata and H. Nunokawa, Phys. Lett. B495 (2000) 369; Nucl. Instrum. Meth. 



A472 (2001) 421; J. Sato, Nucl. Instrum. Meth. A472 (2001) 434; B. Richter, [hep- 
ph/0008222j . 



[23] C. Albright et al, |hep-ex/0008064| ; V. Barger, S. Geer, R. Raja and K. Whisnant, Phys. 
Rev. D 63 (2001) 113011; J. Pinney and O. Yasuda, Phys. Rev. D 64 (2001) 093008. 



[24] Y. Itow et ai, The JHF-Kamioka neutrino project, |hep-ex/0106019 



[25] The MINOS Collaboration, P. Adamson et ai, MINOS Detectors Technical Design 
Report, Version 1.0, NuMI-L-337, October 1998. 

[26] OPERA Collaboration, M. Guler et ai, OPERA: An Appearance Experiment to Search 

for Nu/Mu < > Nu/Tau Oscillations in the CNGS Beam. Experimental Proposal, 

CERN-SPSC-2000-028, CERN-SPSC-P-318, LNGS-P25-00, Jul 2000. 

[27] H. Murayama and A. Pierce, Phys. Rev. D65 (2002) 013012. 

See also, A. de Gouvea and C. Pena-Garay, Phys. Rev. D64 (2001) 113011; V. Barger, 
D. Marfatia, and B. P. Wood, Phys. Lett. B498 (2001) 53; R. Barbieri and A. Strumia, 
JHEP 0012 (2000) 016. 

[28] J. J. Gomez-Cadenas et al. (CERN working group on Super Beams Collaboration), 



hep-ph/0 105297 . 



[29] |http: / /www-off-axis.fnal.gov/loi/loi.ps 



35 



